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Efficient metliods for cliaracterizing tlie performance of quantum measurements are important in 
the experimental quantum sciences. Ideally, one requires both a physically relevant distinguishability 
measure between measurement operations and a well-defined experimental procedure for estimating 
the distinguishability measure. Here, we propose the average measurement fidelity and error between 
quantum measurements as distinguishability measures. We provide protocols for obtaining bounds 
on these quantities that are both estimable using experimentally accessible quantities and scalable 
in the size of the quantum system. This gives an exponential reduction in the resources required for 
complete measurement tomography. We explain why the bounds should be valid in large generality 
and illustrate the method via numerical examples. 



I. INTRODUCTION 

Measurement plays a fundamental role in the quantum 
sciences as it allows one to evolve and extract informa- 
tion from a quantum system. From a practical perspec- 
tive, measurements are imperative for the success of a 
wide variety of tasks from quantum information process- 
ing (QIP) and communication [ ] to performing high- 
precision metrology [ ]. If large-scale QIP is achieved 
experimentally, certain tasks that are conjectured to be 
classically inefhcient will become realizable using a scal- 
able number of resources 4]. 

Various models of computation have been proposed for 
QIP, such as the standard circuit model [ ], measurement 
based quantum computation (MBQC) [G], and topologi- 
cal quantum computation [ ]. Each of these models re- 
lies greatly on the ability to perform accurate measure- 
ments. A large number of physical systems have been 
proposed as candidates for implementing QIP. A short, 
and certainly non-exhaustive, list of systems includes su- 
perconducting circuits [^', Nitrogen- vacancy (NV) cen- 
ters [id], trapped ions [j i], NMR [, quantum dots [i 
and optical implementations [14]. Measurement schemes 
vary greatly across these systems and so it is important 
to have straightforward protocols for comparing measure- 
ments that are both independent of the particular type 
of implementation and scalable in the size of the system. 
Providing such protocols is the main goal of this paper. 

There has been a significant amount of research to- 
wards completely characterizing the error on the opera- 
tions used to process and measure quantum information 
via quantum process [15-17] and measurement tomog- 
raphy [IS, 1!)]. In principle, these tomographic meth- 
ods can be used to characterize the error affecting any 
quantum operation or measurement realized in an exper- 
imental setting. Unfortunately, there are various signif- 
icant drawbacks to complete process and measurement 
tomography. Process tomography requires an exponen- 
tial number of resources in the number n of quantum 
bits (qubits) that comprise the system (the number of 
parameters required to just describe the process scales 



as 16"), is not robust against state-preparation and mea- 
surement errors [2(1], and requires intensive classical post- 
processing of measurement data [21]. Not surprisingly, 
since measurement tomography is realized by essentially 
inverting process tomography [18, 19], it suffers from sim- 
ilar drawbacks to quantum process tomography. In par- 
ticular, complete measurement tomography requires the 
ability to prepare a complicated set of pure input states 
with extremely high precision (possibly by performing 
complex unitary gates), scales badly in the size of the 
system, and can require lengthy post-processing of the 
tomographic data. 

In many cases, one may only be interested in a sub- 
set of parameters characterizing the noise, or determin- 
ing the strength of the noise process rather than the en- 
tire process itself. For instance, in fault-tolerance [22- 
25], as long as the strength of the noise affecting state- 
preparation, gates, and measurements is bounded by 
some value, large-scale computation is possible. As a re- 
sult, various methods for partially characterizing the er- 
ror of quantum gates [20-35] have recently been proposed 
which attempt to circumvent many of the problems asso- 
ciated with full process tomography. The error is defined 
via the "quantum average gate fidelity" , which is derived 
from the quantum channel fidelity between two quantum 
processes. Since quantum gates are unitary operations, 
the average gate fidelity (and higher order moments) take 
a simple form and can be calculated analytically [36-38] . 

Similarly, partial estimation of measurement errors by 
more efficient methods than full measurement tomogra- 
phy would be extremely valuable, but have yet be consid- 
ered in significant detail. There are a variety of questions 
to consider in this endeavor, such as how to deal with 
abundance of different measurement schemes across the 
various implementations, and the ambiguity regarding 
what the single parameter characterizing the error should 
be. Since measurements are non-unitary processes, the 
quantum channel fidelity between the measurement and 
noise process is a very complex quantity and cannot be 
calculated analytically as is the case for quantum gates. 

Here, we propose the average measurement fidelity and 



2 



error rate of a quantum measurement as natural mea- 
sures of the error, and provide scalable, experimentally 
implementable protocols for estimating these quantities. 
The output of the protocol is an upper (lower) hound on 
the error rate (average fidelity) of the measurement. We 
focus mainly on the case where we attempt to implement 
a finite-dimensional projective measurement (also known 
as a projective- valued measure (PVM)), but instead im- 
plement a general quantum measurement in the form of a 
positive operator- valued measure (POVM). In particular, 
we analyze the following two scenarios: 

1. The ideal measurement is a rank-1 PVM and we 
are only concerned with characterizing measure- 
ment probabilities. 

2. The ideal measurement is a rank-1 PVM and we 
are concerned with simultaneously characterizing 
measurement probabilities and output states. 

The case of rank-1 PVM's is significant in QIP since 
each computational model mentioned above can achieve 
universal quantum computation using rank-1 PVM's. We 
leave extending our results to arbitrary- rank PVM's as 
a direction of future research, however we anticipate the 
results presented here can be carried over to these mea- 
surements. Moreover, while the most general measure- 
ment model is a POVM, any POVM can be realized as 
a PVM on an extended Hilbert space by attaching an- 
cilla systems to the Hilbert space of interest [ ]. Thus, 
characterizing the quality of PVM's can provide direct 
information about one's ability to implement POVM's. 

The methods we present here to obtain bounds on error 
rates of measurements are general and should be straight- 
forward to implement in most systems. The protocol in- 
volves the ability to prepare a random state from the ba- 
sis that constitutes the ideal measurement, and perform 
the noisy measurement (perhaps consecutively). Com- 
puting the upper bound from the measurement data is 
straightforward, and generic estimates of how many tri- 
als need to be performed can be obtained using statistical 
methods. We provide a direct comparison of the time- 
complexity between our protocol and that of completely 
reconstructing the noisy measurement POVM elements 
via tomography. Our method scales independently of d 
(where d is the dimension of the system, e.g. d = 2" 
for an rt-qubit system) while a complete reconstruction 
scales as 0{d?). Thus, while our method gives much less 
information than full tomography, it is scalable, straight- 
forward to implement, and involves less post-processing 
of measurement data. 

The structure of the presentation is as follows. First, 
in Sec. (II) we set notation, discuss the framework of 
PVM's and POVM's, and define the noisy measurement 
models. Next, in Sec. (Ill), we define the quantum average 
measurement fidelity and quantum average measurement 
error, which are the operationally relevant metrics we use 
to compare measurement operations. 

In Sec. (IV) we present the main result for the first of 
the two scenarios listed above. We provide an expression 



for a lower (upper) bound on the fidelity (error rate), 
and present an experimental protocol for obtaining these 
bounds. In Sec. (V) we give a derivation of the bounds 
and show they satisfy certain necessary conditions to be 
useful in practice. Sec. (VI) provides a discussion of the 
general validity of the bounds and numerical examples 
for the case of a single-qubit system are analyzed. Our 
derived bounds are shown to be valid for every example 
analyzed, confirming the discussion regarding the validity 
of the method. 

We discuss the extension of the method to rank-1 
PVM's with output states in Sec. (VII) and provide the 
experimental protocol for this case in Sec. (VIII). The re- 
source analysis for the protocols are provided in Sec. (IX) 
and concluding remarks are made in Sec. (X). The reader 
interested in the bounds and protocols without most of 
the technical details is referred to the text surrounding 
Eq. (4.2) for rank-1 PVM's and the text surrounding 
Eq. (7.14) for rank-1 PVM's with output states. 

II. MEASUREMENT OF QUANTUM SYSTEMS 
AND THE ERROR MODEL 

Denote the quantum system by S and suppose it is 
represented by a Hilbert space % of dimension d < oo 
[d — 2" for an n-qubit system). The most general type 
of measurement one can make is a "positive operator- 
valued measure" (POVM). A POVM consists of a set 
{Ek} of linear operators on % that satisfy, 

Ek > 0, 

Y^Ek = l. (2.1) 

k 

If the input state to the measurement is a, the probability 
of obtaining outcome "fc" is given by 

rk = tr {Eka) . (2.2) 

An important subset of POVM measurements are 
"projection- valued measures" (PVM's), which corre- 
spond to the case of each Ei^ being equal to a projection 
operator H^. Hence, in addition to the conditions listed 
in Eq. (2.1), H^ = Hfc for each k. 

PVM's are a very general class of quantum measure- 
ments and, as previously mentioned, are sufficient for 
performing universal quantum computation in various 
computational models [5, G]. For instance, computa- 
tional basis measurements allow for universality in the 
standard circuit model and single-qubit projective mea- 
surements provide universality in MBQC. Other impor- 
tant examples of PVM's in quantum information theory 
are parity measurements, which are used extensively in 
quantum error-correction and fault tolerance [40-43]. It 
is important to note that, by Naimark's theorem [39], 
any POVM can be performed via implementing a PVM 
on a larger Hilbert space. Hence, methods for determin- 
ing error rates of PVM's can also give direct information 
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regarding the error associated with performing general 
POVM measurements. 

PVM's are in a 1-1 correspondence with observables 
(Hermitian operators) O, where O is non-degenerate if 
and only if the PVM consists of rank-1 projective ele- 
ments . We write the observable O as 



(2.3) 



fc=i 



where the eigenvalues correspond to the measurement 
values and 11^ := X]j=i is the rank-dfc projec- 

tor onto the eigenspace spanned by the eigenvectors 

Hence, X]fc=i dk = d and the PVM is rank-1 if and only 
if dfc = 1 for every k. The ideal quantum operation 
that corresponds to the PVM (without post-selection), 
denoted A4, has action on input state cr, 



X(a) ^^n^aHfe =^Pfc 



fe=i 



fe=i 



tr (HfefT) 



where we define, 



Pk = Pkicr) = tr (life cr) . 



(2.4) 



(2.5) 



In the case of analyzing only measurement probabil- 
ities, the PVM is a mapping from the set of quantum 
states to the probability vectors {pi, ...,pi,), 



Mia) = (pi(a),...,pb(a)) = {tr (Uia) ^ 



..,tr (Hfeo-)). 

(2.6) 



This case is especially interesting since in many scenar- 
ios one is mainly interested in the output value of the 
measurement rather than the output state itself. For 
instance, the output of a quantum algorithm is usually 
the value of a measurement in the computational basis. 
Hence, as long as the measurement values are obtained 
with correct probabilities, the measurement is deemed 
successful. 

In addition, in many current implementations of mea- 
surements, the action of performing the measurement op- 
eration to obtain the output value destroys or drastically 
alters the state of the system. For instance in photo- 
detection, which forms the basis of various measurement 
schemes in optical, atomic, and superconducting systems, 
the measurement consists of recording the occurrence of 
a photon. If the state of the system is encoded into a 
degree of freedom (mode) of the photon, such as polar- 
ization or frequency, detection of the photon can record 
the output of a measurement of this degree of freedom. In 
common photodetectors, such as avalanche photodiodes, 
the photon is lost in the process of creating a current 
through the photoelectric effect. Hence, the system en- 
coding the information to be measured is destroyed, but 
the measurement output can be accessed. 

An example where the system is not destroyed yet 
the state is not preserved from the measurement is the 



Nitrogen- Vacancy (NV) center in diamond [10]. In this 
case, the processes of measurement and ground-state po- 
larization are identical. Hence, while one can obtain the 
relevant measurement outcomes and statistics, the out- 
put state of the measurement is always the ground state. 

In the general case of Eq. (2.4), the noisy measurement 
£ is modeled by 



b 

E 

k=l 



rk{a)Pk{o-) 



(2.7) 



where we allow both the noisy measurement probabilities 
rfc and output states pk to be functions of cr. By Eq. (2.4), 
rfc(cr) and /Ofc(cr) are ideally given by Pfe(cr) = tr(nfccr) 
and t^^n^*) respectively. If we analyze only measurement 
probabilities, we have 

E{(t) = (ri(a), ...,n{a)) = (tr {E^a) , tr {E^a)) (2.8) 

where £ is allowed to be of a completely general form 
by assuming it is modeled by a POVM {E^^}\^^. The 
first scenario we are interested in is to compare the prob- 
ability distributions in Eq. (2.G) and (2.8). We analyze 
the more general case of measurement probabilities and 
output states (ie. comparing Eq.'s (2.4) and (2.7)) in 
Sec. (VII). We first discuss the figure of merits we will 
use to compare the ideal measurement process M and 
noisy measurement process £ . 



III. QUANTUM AVERAGE MEASUREMENT 
FIDELITY AND ERROR 

A completely natural question is, how should we com- 
pare the ideal and actual measurements and £1 A 
set of criteria that a distance measure. A, for comparing 
ideal and real quantum processes should satisfy has been 
given previously [44]. In practice, no A will satisfy all 
of these criteria simultaneously. Thus, one must settle 
for A to satisfy a subset of these criteria, in addition to 
other criteria that may be useful for the particular task 
at hand. 

As mentioned previously, the average gate fidelity J-gfA 
is a useful method for comparing an intended unitary 
operation hi and actual quantum process £. There are 
various reasons for the utility of the average gate fidelity, 
for instance, it satisfies the following properties: 

1. There is a straightforward method for evaluating 
J-£ u (given a description of hi and £), 



2. J-£_u h&s a well-motivated physical interpretation, 

3. All states are taken into account in an unbiased 
manner when calculating J-gu^ 



4. J-£pi is experimentally accessible via efficient pro- 
tocols. 
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An important drawback of the average gate fidelity is 
that it is not a metric. We would like similar properties 
to hold for our method of comparing ideal and real mea- 
surements. Let us briefly outline how the average gate 
fidelity is derived from more general quantities, which 
will provide intuition for how to define our method for 
comparing measurements. 

The average gate fidelity is derived from the state- 
dependent quantum channel fidelity, which is a standard 
method for comparing quantum operations. If £i and £2 
are quantum operations and a is a quantum state, the 
quantum channel fidelity between £1 and £2, denoted 
Fs^^S^, is given by the standard state fidelity between 
£i{i7) and £2(0"), 

Fs,,£,{a)^F{£,{cy),£2{a)) 

\v^^/£^(a)£^{&)^/£^)^ . (3.1) 

When one of the operations in Eq. (3.1) is unitary (say 
£2 =U), the channel fidelity is called the quantum gate 
fidelity, and when a is pure (a = \ip){'ip\)^ the gate fidelity 
takes the extremely simple form, 

Fe,,umm = ^^{£immum{m 

= {i^\m)mi^) (3.2) 

where 

k=U^o£. (3.3) 



The average quantum gate fidelity, denoted J-g-^ u j is ob- 
tained by integrating over all pure input states. The 
integral is taken over the unitarily invariant Haar mea- 
sure (also known as the Fubini-Study measure) on the set 
of pure states [-! "']. In this paper we denote the Fubini- 
Study measure by /i. This gives. 



(^|A(|7/^)(V'|)|^)dV 
E,tr(A,)tr(A,-) + d 



(3.4) 



where {Aj} is any set of Kraus operators for A [ : ]. Thus, 
the average gate fidelity reduces to an extremely simple 
form because one of the operations is unitary. 

Following this intuition, from Eq.'s (2.4), (2.7), 
and (3.1), we have that for the ideal (A^) and noisy {£) 
measurements , 

F{£{a),M{a)) 




where the state-dependence in pk, rj-, pk is omitted 
for notational convenience. Unlike the unitary case 
(Eq. (3.2)), such an expression does not reduce to a 
simple form in general. However, ideally, the figure of 
merit we use to distinguish M and £ will contain evenly 
weighted information from all possible input states and 
also have direct operational significance. As a result, 
we define the quantum average measurement fidelity, de- 
noted Fg^M, to be the Haar integral of Eq. (3.5) over 
pure input states. 



F£,M = J Fi£m),Mm))- (3.6) 
We also define, 

r£,M = 1 - Fs^M (3.7) 

to be the quantum average measurement error. In the 
case of only being concerned with measurement proba- 
bihties, that is M and £ are given by Eq.'s (2.6) and (2.8) 
respectively, we have by Eq. (3.5), 



F£,M = J VPkrkj dip 



Kk=l 



(3.8) 



rs^M = 1 - / I E Vt^(nM¥\MEMW)j d^. 

(3.9) 



The direct relationship between rg j^ and Fg im- 
plies statements and bounds proven about one directly 
applies to the other. Since Fg^M and r^jvi will not have 
simple forms like the average gate fidelity in Eq. (3.4), our 
goal is to provide efficient methods for estimating rg^M 
(equivalently Fg^j^) that can be experimentally imple- 
mented in a simple manner. 

An important reason for using the average measure- 
ment fidelity (error) as defined above to distinguish M 
and £ is that this provides a state-independent distance 
measure to compare the operations which can be con- 
nected to the diamond norm distance [46] between quan- 
tum operations [47]. Computing the diamond norm is 
an exponentialy hard task since one needs a complete 
description of the quantum operations. Thus, the dia- 
mond norm is neither straightforward to calculate nor 
experimentally accessible. However, it is commonly used 
in fault-tolerant analyses for determining threshold error 
rates of physical operations. Thus, information about 
the diamond norm provided by the average measurement 
error defined above can potentially provide information 
regarding fault-tolerant computation. 
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IV. EXPERIMENTAL PROTOCOL: RANK-1 
PVM'S 

In this section, we look at obtaining an upper bound 
on r£^M in the case of rank-1 PVM's where we are not 
concerned with output states. Again, we note that upper 
bounds on jvi are equivalent to lower bounds on Fg^M 
and so for clarity we phrase the following discussion in 
terms of r^,^. 

There are d PVM elements, {Hi, Xld}, each of which 
is a rank-1 projection operator. From Eq. (2.8), the out- 
come of the measurement of a state can be asso- 
ciated to a member from {1, d} with frequency distri- 
bution 

(ri,...,rd) = {iv{Ei\i;){i^\),...MEd\^){m- (4.1) 



Goal : Obtain an upper bound, uh, for ts^m as defined 
in Eq. (3.9). 

The experimental protocol to obtain ub is as follows: 
Protocol : 

Step 1: Choose a pair of indices {l,m) uniformly at ran- 
dom from the set V = {0, d — 1} x {0, d — 1}. 

Step 2: For each j g {/,m}, 

a) : Prepare the quantum state perform the noisy 
measurement £ on Hj = and record whether out- 
come "j" is obtained, 

b) : Repeat a) many times and denote the frequency of 
obtaining by Uj, that is, Uj — tr{IljEj), 

(see Sec. IX A for a discussion of the number of repeti- 
tions required to estimate Uj to a desired accuracy and 
confidence) . 

Step 3: Repeat Steps 1 and 2 K times, where K is dic- 
tated by the desired accuracy and confidence in estimat- 
ing ub 

(see Sec. IX B for a discussion of the size of K). 

Step 4: Compute the upper bound of rg^M, ub, via the 
formula (derived in Sec. V), 



ub^l- 



1 + d 



1 + dX 

1 + d 

{1-X) 



where 
X 



(4.2) 



(4.3) 



and {(li, I2), (-K'l, -ft'2)} are the K trials dictated by 
Step 3. 

This concludes the protocol. 

There are various important points about the protocol 
that should be emphasized. First, the number of trials 
required in Steps 2b) and 3 are independent of d. Thus, 
the time-complexity of the entire protocol is independent 
of d, and depends only on the desired accuracy and confi- 
dence of the estimate of ub (see Sec. IX). Second, ub can 
be computed from the above protocol using only: 

1. Applications of the noisy measurement and, 

2. The ability to prepare the d pure input states \ipj). 

Lastly, it is straightforward to show the following two 
properties of ub, (see Sec. VD) 



1. In the limit of r£,M i 0, 



ub \, 0, 



(4.4) 



2. ub scales well in d. 



These are clearly necessary conditions for ub to be a good 
upper bound on r^^^vi. Moreover, property 1 implies that, 
loosely speaking, in the small error limit ub can be taken 
as an estimate of r^ jvi. 



V. DERIVATION OF THE UPPER BOUND 

From Eq.'s (3.7), (4.2), and (4.3), showing ub is an 
upper bound for r^ j^ is equivalent to showing 



lb 



did+1) 



(5.1) 



is a lower bound for Fg^M- First, we have from Eq. (3.8), 
F{£{<j),M{<7))^(^^F^^ 

(5.2) 

where the state-dependence is implicit in the pk and ■ 
Taking the integral over all pure states gives. 



F{£{\^){^\),M{m{md^ 



E 

fc=i 



rkPkdtp 



E 



^/rirmPiPnidip 



(5.3) 



To proceed, by Schur's Lemma, we have that for any 
positive integer t, 



{kiM) 



\iP){i;\'^'dtlj 



Ilsyniit,d) 

tr [nsyin(i, d)] 



(5.4) 
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where Hgyrait, d) is the projector onto the symmetric sub- 
space of the t-partite Hilbert space whose factor spaces 
are each c?-dimensional Hilbert spaces [38, 48]. It is 
straightforward to show that Ilsyni{t,d) is equal to the 
normalized sum of the t! elements in the group of permu- 
tation operators on the t-partite Hilbert space. We look 
at the sums in Eq. (5.3) separately. 



A. X] / '''kPhdij} 

We have from Eq. (5.4) 

rkPkdi}j= I iT{Ek\i^){-^\)ir{IVk\-^){i}j\)d^ 



ti{[Ek®nk]\4j){^\®\ij){'4,\)d',p 

2ir{[Eu®Iik\^syu.{2,d)) 



d{d+l) 



(5.5) 



Since the group of permutations of two objects consists 
of the identity and SWAP operations, we have 



nsym(2,d) 



1(E}1 + SWAP 



(5.6) 



Therefore, 



rkPkdip 



tr(SfeHfe) + dtr(i;fe|) 



d{d+l) 



(5.7) 



and so, since X]fc=i -^fc = 1; 



^ [ ^, ^ tr(^fcnfe) + dtr(£;4) 
2^ / rkPkdi^ = 2^ - 



k=l 



k = l 



d{d+l) 



T,k=l "fc + d 
d{d+l) 



(5.8) 



which gives the first two terms in the numerator in 
Eq. (5.1). 



B. ^ J y/rirmPiPmdi> 



We have, 



\/rirmPWmdil) 

Vtr(n,|V')(^|)tr(S,|V;)(V|)Vtr(n,„|V^)(V;|)tr(i;„.|V)(V'l)rf^. (5.9) 

Computing the above integral is analytically is not possible because of the square root in the argument. Now we 
assume the following inequality holds (see Sec. VI), 



V'tr(H,|V')(V'|)tr(i?j|^)(V'|)\/tr(H„|V)(VI)tr(i;™|V)(^|)dV' 
> / v/tr(n,|^)(V^|)tr(n,i;jH,(|V)(^Al))v/tr(n„,|^)(7^|)tr(H,„£;„H„(|V')(V'|))dV 
= j V(V'/|^^/|V'i)tr(H,(|V')(V'|))V(V'm|^^m|V'm)tr(H„,(|i^)(V'|))rfV 

= / tr(H,(|V')(VI))tr(H„(|V')(V'|))dV. (5-10) 



Note that for each j it is generally not true that Ej — ujUj From Eq. (5.4) we have 

is positive semidefinite. If this were the case then the 
above inequality would always hold. 



tr(H;(|^)(V|))tr(H„(|7^)(V|))dV' 



= / tr[{Ili(g>n„,)\ib){^\(g>\,p){t^\]di^ 

tr (n;H„)+tr (HQtr (H^) 
d{d+l) 
1 

~ d{d + l) 



(5.11) 
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since 11/ and are projectors onto orthogonal sub- 
spaces. Thus, 



Vtr(^;n;)tr (^„n„) 
d{d+l) 



y/UlUm 
d{d+l) 



(5.12) 



which gives the last term in the numerator of Eq. (5.1). 
In total, assuming the inequality in Eq. (5.10) holds, 



1 + dX _d + J2k=i "fe + y/m 



1 + d 



d{d+l) 



(5.13) 



is a lower bound for Fg^M- Before analyzing the validity 
of the lower bound, we briefly compare its expression 
with that of the average gate fidelity and show that it 
satisfies certain necessary conditions for it to be useful. 



C. Comparison With Quantum Gate Fidelity 

There is a nice parallel between the lower bound on 
the average measurement fidelity and the exact expres- 
sion for the average quantum gate fidelity. As previously 
mentioned, the quantum average gate fidelity takes a par- 
ticularly nice form because one of the operations is uni- 
tary [■'' ]. Indeed, if one compares the unitary U and 
quantum operation £, the average gate fidelity between 



U and E, J-£.u, is given by, 



£.U 



Ad+l 
d+l 



where 



(5.14) 



(5.15) 



and {Ak\ is any set of Kraus operators for the quantum 
operation h.=U^ o £. 

Our expression for a lower bound on the average mea- 
surement fidelity Fg^M takes a similar form. 



Ih- 



Xd+l 
d+l 



where 



x = ^Yl v/tr(ni^Otr(n,„£;,„). 



(5.16) 



(5.17) 



In some sense this is not surprising since each quantity is 
a Haar integral over functions of two copies of a quantum 
state \ip){tp\. A direction of further research is to under- 
stand properties of the average measurement fidelity in 
more detail, and draw more parallels with well-known 
distinguishability measures such as the average gate fi- 



delity. 



D. Necessary Conditions For ub To Be a Useful 
Upper Bound on r£,7\4 

In this subsection, we discuss two necessary conditions 
lb must satisfy in order to be a useful lower bound on 
Fg.M- By the equivalence between ub being an upper 
bound for rgj^ and lb being a lower bound for Fg j^, 
these conditions are equivalent to corresponding condi- 
tions for ub. 

Two necessary conditions that lb must satisfy in order 
to be a useful lower bound on Fg ji^ are, 

1. Limit of No Error: As the measurement error goes 
to 0, lb t 0, and, 

2. Scaling in d: lb scales well in the dimension d of the 
system. 

Here, we show lb satisfies both of these criteria. 



1. Limit of No Error 

Suppose the measurement error goes to in that the 
POVM elements converge to the ideal PVM elements, 



This implies for each j e {1, ...,cJ}, 

Uj = tr(7rj£'j) t 1- 



(5.18) 



(5.19) 



Hence, from the set of equations in Eq. (C.l), for each 
{l,m), 



which implies. 



Ibtl- 



(5.20) 



(5.21) 



Thus lb satisfies the necessary condition of converging to 
1 in the limit of no errors. 



2. Scaling in d 

Let us now show that lb scales well in d. By Eq. (6.2), 
suppose that for each {l,m), ei^m > is such that. 



(5.22) 
This is equivalent to, 

(5.23) 



d{d + l) 
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Since e;.„i = if Z = to, summing over all I m gives, 

5: / f,Mm - E + E (5-24) 



Assume the error in the POVM elements remains roughly 
independent of d in that tr{IljEj) is independent of d. 
Let us suppose for the sake of clarity of the argument 
that 

trinjEj) = 1- S (5.25) 
for some small 6 > independent of j and d. Then, 



E 



^ d{d+l) 



{l~6)d{d-l) 
d{d+l) 



(5.26) 



Hence as d — >■ cx) 



Define 



E / ^ (1 - 5) + E (5.27) 



e = E ^'.^ 



Then, 



(5.28) 



(5.29) 



and so, since F^^j^ < 1, each value of ti^m must scale 
as Hence, Z6 scales well in large dimensions since it 
essentially depends on the size of the Uk (ie. 1 — 5). 



VI. VALIDITY OF THE UPPER BOUND 

The upper bound in Eq. (4.2) is valid provided the in- 
equality in Eq. (5.10) holds. The goal of this section is to 
show that this inequality holds in very general situations. 
To set notation, we define 



fiM) = v/tr(Hz|V')(V'|)tr(iJ,|V)(V'l)\/tr(H™|V^)(V'|)tr(i5„|V)(^|), 

9iA\i^)) = v/tr(Hz|V')(V'l)tr(H;£;ini|V)(VI)v/tr(n„|V')(V|)tr(n„£;„H„|V)(V'|) 
= V7/Iii;;:tr(H,|V')(V'|)tr(H™|V')(^Al) 



(6.1) 



so that we want to show, 

/ [kmm)~9l,m{\mdi^>Q- (6.2) 

We have, 

[kM) - 9i,M)]d'>P ^ j ki,nm)hi,M)di^ 

(6.3) 

where 



or 



ki,m{m v/tr(ndV')(V'l)tr(n™|^)(V;|) (6.4) 

KM) ■■= Vi^{EM¥¥^iKMW) (6.5) 
- y^II^v/tr(H,|V')(V'l)tr(H,„|V)(V'|). (6.6) 

Since fc;_m(|?/')) > for all we have, 

kmm) - 9Lmi\^)) < ^ KM)) < 0. (6.7) 

However, /ii,m(|V')) < if and only if one of 

{MEi-uiUim<0 (6.8) 



\E„ 



1. 



< 



(6.9) 



is true. Without loss of generality, suppose {ip\Ei — 
w/H/|?A) < holds. 

As a matrix written in the {|^fe)} basis, we have 
El — uiHi is equal to Ei except in the {I, Z)'th entry which 
is 0. If IV') = \^i) then {i}\Ei - mzH/IV-) = 0, and if is 
orthogonal to \ipi), {iplEi-uinilip) > 0. Hence, Eq. (6.8) 
can only be satisfied if {ip) has large enough amplitude in 
the subspace defined by {ipi) . The amount of amplitude 
that is required depends on the size of the off-diagonal el- 
ements (coherence) in the Vth column (row) of Ei. If the 
coherence is not too large then \ip) will need to have large 
amplitude in \i{ji), and as the coherence grows larger, \^p) 
can potentially satisfy Eq. (6.8) without being close to 

The key point is that, if \^) has large amplitude in 
IV';), and so there is a chance of {ip\Ei — uiUilip) < 
0, then tr(Hm|V')('/'l) is small. Thus, fc;, ,„(]?/')) will 
be small and so fcj^„i(|V')) acts like a modulating fac- 
tor in the expression for fi.m{\ip)) — 5i,m(|V')) to ensure 
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fi,m{W) - 9i,m{\'^)) is small when lip) is close to \tjji) (or 
\4'm) by symmetry). Moreover, the measure of the set 
of states for which \^p) is close to lipi) (or |?/'m)) is small, 
and by Levy's Lemma [-1')] decreases exponentially in the 
dimension d. Thus, we have that the measure of states 
for which \ip) is close to (or \tpm)) is small and, for 
such a state, fi,mi\'4')) ~ 9i,mi\''P)) is small in magnitude. 

In total, if the POVM element Ei does not have strong 
coherence then, since {ip) must be close to \4>i) to make 
fi,m - gi.ni < 0, we expect 



Since our goal is obtaining a lower bound on J (/o,i — 
go^i)(\ip))dij^, let us attempt to understand when one of 
Eq.'s (6.8) or (6.9) is true (ie. {iplEo - UoUo\ip) < or 
{ip\Ei — uilli\ip) < 0) since this is a necessary condi- 
tion for hi^,n{\^p)) < 0. Without loss of generality, let us 
assume 



\Eo - uolloli^) <0. 



(6.15) 



We have, 



'Amim - 9Lmm)]d^ > 0. (6.10) ('^l^o - noUom = 2Re (^7 cos (J^ e'^sin (^0 



Hence, our bounds should be valid for most practical 
cases since, realistically, the coherence will not be over- 
whelmingly large. Note that the bounds are guaranteed 
to be valid if there is no coherence at all, that is, if the 
POVM elements are diagonal in the {[V'fc)} basis. In the 
next section we perform the above analysis for a single 
qubit and in Sec. VI B we numerically investigate the 
single-qubit case. Interestingly, we find that for all val- 
ues of coherence magnitude, the bounds are valid. This 
provides evidence that the bounds will be valid in large 
generality. Clearly, a more in-depth investigation of suf- 
ficient conditions for the bounds to be valid is desirable, 
however the above argument and numerical results of 
Sec. (VI B) indicate that the bounds and protocol will 
be valid in most practical cases. 



A. Single-Qubit 

Our goal in this section is to obtain a lower bound for 



(6.11) 



and provide sufficient conditions for when Eq. (6.2) holds 
in the case of a single qubit. We use the Bloch sphere 
representation of a single qubit: 



cos 



e"^ sin 



ID 



(6.12) 



where 6 e [0,7r), e [0,27r). As the results will depend 
on the magnitude of the coherence in the POVM ele- 
ments, we assume without loss of generality that 7 G R 
and 7 > 0. 



7 



7 tY{Eani 



^ ii[EiU„) -7 

1 I -7 Ml 



(6.13) 



(6.14) 



-Ftr(£;oni)sin^ 

and it is straightforward to show ("01 -Eq - 
satisfied if 

vr Stt 



(6.16) 



mqIIoIV') < is 



2' 2 



0, 2 arccot — 



tr(i^IIi) 
27 cos((/)) 



(6.17) 



Hence, denoting the set of states that satisfy ('0|£'O — 
moHqIV') < by Ao,7, we have that Ao,-y is contained in 
the set of all states with and 9 given by Eq. (6.17). 
As expected, the measure of ^0,7; A'(^o,7): is extremely 
small for weak coherence and grows larger as the coher- 
ence increases in magnitude. 

One can obtain an upper bound for /i(v4o.^) by noting 
that the uniform measure on the Bloch sphere has density 
function 



1 

47r 

Hence, /i(Ao,-y) is at most 



sin(6l). 



1 
4 



1 

47r 



cos 2 arccot 



tr(£;oHi 



27 C0S((/)) 

Now, if IV) e Ao,^, e < 2 arccot 

tr(£;oHi) 



(6.18) 



d(j). (6.19) 



and so, 



tr(Hi|'(/')(?/'|) < sin^ arccot 



27 cos((/)) 



COS^((/)) 



cos2(0)+(*^ 



(6.20) 



Hence, we have that for weak coherence 7, the set of 
states which satisfy Eq. (6.15) will have small measure 
(bounded by Eq. (6.19)) and by Eq. (6.20), the value of 
5o,i on these states is small. Using a symmetric argument 
for the state |1) gives the following set of equations. 
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^(^1.7) < 

If |V)GAo,7:tr(ni|V')(VI)< 



1 1 

4 47r 

1 1 

4 47r 



/''cosf2arccotf-^ 
if V V 27Cos(<^ 



cos 2arccot 



cos^((?!)) 



If |V^)eAi,^,tr(no|^)(^|)< 



cos^((?!)) 



with e 



tr(£;ino) 

27 C0S((/)) 

2' Y 



COS 



with d) G 



2' T 



(6.21) 



Using these results we can obtain a lower bound for Hence, in total, 
/Vi(l^))rfV' = /(/o,i-3oa)(l^>W- We have. 



(/o,i - go,i){m)d^ - E / (/o-i - .9o,i)(lV'))dV' 
(/o,i-5o,i)(l^))dV'. (6.22) 



'(Ao.tUAi,^)^ 

Using the set of equations in Eq. (6.21) one can show 



> 



(/o,i-3o,i)(IV')W 

7 

- v^tr(no|v)(^|)tr(ni|?/.)(^|)dv 



27r 



By Eq. (6.20), 



. 4 / / tr(^ino 
sm arccot — 



27 cos(0) 



d(t). 



> 



2tt 



and similarly. 



(/o,i"3o,i)(|V'))rfV' 



COS'^((/)) 

cos2(0)+(*^ 



(/o.i-3o,i)(|V'))rfV' 



(6.24) 



> 



27r 



COS^((/)) 



cos2((/)) 4 



tr(gino) 

27 



(/o,i-5o4)(l^))c«^ 



> -(<50,7 + (^l^-y) + 



(Ao,tUAi,^)<= 



(/o,i-5o,i)(|V'))d^ 



(6.26) 



where 



/i((Ao^^U Ai,^)'^) 

1 1 



> 



2 47r 
1 

47r 



cos 2arccot — 



cos ^2arccot 



tr(£;oni) 

27 C0S((/)) 

tr(^ino) 
27 COs((/i) 



3.27) 



(6.23) ^"^^ 1^) ^ (^o.-r ^ ^1.7)''' 

[/o,i-.9o,i](|^))>0. (6.28) 
Eq. (6.26) provides an explicit lower bound on 

/ /io,i(|i/'))rfV' = / (/o,i - 9o,i)iH))di}. Next, we analyze 
numerical results for the single-qubit case 



B. Numerical Examples 

We now look at some simple single-qubit examples to 
help clarify the technical details presented thus far. In 
order to perform integration over the Fubini-Study mea- 
sure on single-qubit states, we use the Bloch sphere rep- 
resentation of a single qubit state. 



cos - |0)+sin - e'^ll) 



(6.29) 



(6 25) where 6* G [0,7r], </> e [0, 27r]. Fubini-Study integration 
reduces to uniform integration over the unit sphere C 
R'^. For any function f{d,4>) we have that the integral 
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of / with respect to this uniform probabihty measure is 
the double integral, 

-T / f{0,(t))sm{9)ded<l>. (6.30) 



With Eq. (6.32) in hand, we look at the following exam- 
ples. 

Suppose our goal is to characterize a single-qubit pro- 
jective measurement in the computational basis. The 
POVM elements {Eq, Ei} that model the noise must each 
be positive semidefinite and Eq+Ei = 1. Hence, we write 

where 7 S (D is given by 

7 = Re"'' (6.35) 

and (p e [0, 27r). We have the linear operator Eq (Ei) is 
positive semidefinite if and only if its leading principal 
minors are non-negative. Since d ~ 2, this reduces to the 
condition, 

I7I - i? < v/wo(l-Mo)- (6.36) 

For various cases, we observe how the exact measure- 
ment fidelity compares to uh as R approaches its maxi- 
mum value. Since the results will only depend on R and 
not (/), without loss of generality, one could assume j — R. 
As a verification of this, the numerics given below were 
performed for various G [0, 27r) and, as expected, the 
results were independent of 0. 



The normalization ensures we have a probability mea- 
sure. We have the following quantities. 




Pi = 1 - Po, 

ro = cos^ (2) + ^^0^' sin^ (J^ + Re (iJ^ e'^) sin(0), 
ri = l-ro. (6.31) 
It is straightforward to show. 



..32) 
I 

Three of the cases we analyzed are Uq — 0.99, 0.995, 
and 0.999. The values of lb and ub for each case are 
contained in Table I. Plots of the average measurement 



TABLE I. Values of lb and ub for 3 Cases of uq. 



Uq 


0.99 


0.995 


0.999 


lb 


0.9933 


0.9967 


0.9993 


ub 


0.0067 


0.0033 


0.0007 



error rg^M SiS a function of I7I for each case are given in 
Fig. 1. ' 

Various important features are evident from Fig. 1. 
First, ub is always an upper bound for rg^M (the curves 
in the plot do not intersect). This agrees with the dis- 
cussion in Sec. (V) where we noted that ub should be 
a valid upper bound in large generality (even for large 
coherence) due to both the anti-correlation of the ran- 
dom variables tr(no|V') (V'l) and tr(ni|i/))(i/)|), and the 
measure of the set of states for which 50,1 > fo,i be- 
ing small. Second, as expected, the difference between 
ub and rg^j^ grows smaller as the magnitude of the co- 
herence increases. Thus, the bounds become more tight 
as the coherence in the POVM elements increases. 

These results suggest that the bounds are always valid 
regardless of how large I7I is relative to uq. We tested 
whether there were any violations when uq was varied 
from 0.9 to 0.9999 in increments of lO^'', and the coher- 
ence was varied from its minimum to maximum in each 
case. We found no violations of the upper bound. 



j y/PoroPindip 



1 

An 



J4>=o 



i?°-°COs4 



sm 



E^'^ cos2 



sm 



2Re ( En'^ei^ ] cos 



sm 



1 - E^'^ cos2 



Eq^ sin^ 



2Re{E^'^e^'l>] cos 



sm 



^siii{9)d0. 
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Measurement Error and Coherence Strength 
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0.01 



0.02 



0.03 



Magnitude of Coherence 



a. Note that this assumption is only vaUd because the 
ideal measurement is a rank-1 PVM. We assume the noisy 
measurement process £ is of the form, 



where, as before. 



fc=i 



rfe = tT{Ek(7) 



(7.3) 



(7.4) 



is modeled via a POVM {-E/tj^^i- We assume the pk are 
independent of a since, ideally, the output states 11^ are 
independent of cr. The average measurement fidelity is 
still given by 



F ( ^rkPk,^Pk^k 1 

\ k k J 



dilj (7.5) 



and, as before. 



rsM = 1- / Fs^Mimi^Ddi^- 



(7.6) 



FIG. 1. Plot of exact average measurement error re,M (blue- 
dashed) against coherence strength I7I for three different val- 
ues of uo: 0.99 (upper left), 0.995 (bottom left), and 0.999 
(bottom right). The upper bound ub for re,M in each case is 
also given (green). 



VII. CHARACTERIZING MEASUREMENT 
PROBABILITIES AND OUTPUT STATES FOR 
RANK-1 PVM'S 

Ideally, a quantum measurement produces an output 
state that can be utilized for further purposes. In these 
cases one would like to characterize the error on both the 
measurement statistics and output states of the measure- 
ment. In this section we briefly outline such a method 
for rank-1 PVM's and discuss conditions for the upper 
bound on the measurement error r£ ji^4 to be valid. 

From Eq. (2.4) we have that A4 is given by. 



Mia) = J2^kCTUk = J2Pk^'^ 



(7.1) 



fc=i 



k=l 



where 



Pk = tr (Hfea) (7.2) 
and the output state 11^ is independent of the input state 



In this case where we are also interested in output 
states, the expression for Fg^j^ is much more compli- 
cated. However, we can use strong concavity of the fi- 
delity [511] to obtain. 



FsM = J F{£{\i;)),Mm)) 

> J V^Vtr (pfcHfe)^ d^p. (7.7) 

Since tr (pkHk) is just the fidelity between the fc'th ideal 
and noisy output states of the measurement, we define. 



Fk ■.^Fipk,nk)^tT{pkiik). 



Moreover, since there is only state-dependence in pk and 
rfc, we have. 



F> 



> 



E 



(Lm)e'D 

From Sec. V this gives 



\/nr,nPiP7ndip 



l + dY 



^F^r 



(7.9) 



F£,M > 



1 + rf 



where 



^ E FiF„ 



(7.10) 



(7.11) 
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At this point, we could define lb = (and thus 

ub — 1 — '^i^J ), except the Fk cannot be obtained di- 
rectly from the noisy measurement £. As our goal is an 
experimentally implementable protocol to obtain bounds 
on Fs.M and r^^^, we need to find a lower bound for 
this quantity in terms of accessible quantities. Note that 
the values tr {pkEk) are accessible in that, if the output 
of the measurement is "fc" and one feeds the output state 
Pk into a subsequent measurement, then the probability 
of obtaining output "fc" in the second measurement is 
tr (pkEk). Let us denote Qk to be this probability, 

Qk — IP [obtaining outcome k from £ if the input is pk] 
= tripkEk) (7.12) 

where the notation "P(-)" means "probability of ' . Thus, 
if it is true that for each fc, 

Fk - tr (pfeHfc) >Qk^ tr (pkEk) , (7.13) 

then wc have, 

reM<Y^i^-Z), (7.14) 

where 

X! V";"™ VQiQm- (7.15) 

(Z,m)GX' 

Motivation for why the assumption in Eq. (7.13) will 
often be true is that ideally, for each fc, Ek — Ilk- Hence, 
tr (pkHk) is just the fc'th diagonal element of pk, p\'^ ■ 
Since Ek <l,\i Ek is diagonal in the {|V'fc)} basis, then 
Eq. (7.13) is guaranteed to hold. This is because the 
probability amplitude of Ilk is spread across the diagonal 
elements of Ek- Even when there is coherence in the Ek 
(non-zero off-diagonal elements) we expect the inequality 
should still hold in realistic situations, however a more 
detailed analysis of conditions for this to occur are re- 
quired. Note also that, realistically, one does not need 
Eq. (7.13) to hold for every k, just a large enough num- 
ber to ensure Y > Z . We now provide the experimental 
protocol for determining ub. 



VIII. EXPERIMENTAL PROTOCOL: RANK-1 
PVM'S WITH OUTPUT STATES 

Goal : Obtain ub in the case of imperfect measurement 
probabilities and output states. 

Protocol : 

Step 1: Choose a pair of indices (Z,to) uniformly at ran- 
dom from V = {0, 1} x {0, 1}. 



Step 2: For each j € {I, m}, 

a) Prepare the quantum state Ilj, perform the noisy 
measurement £ on IIj, and record whether outcome "j" 
is obtained, 

b) If outcome "j" is obtained in a), repeat the mea- 
surement on the current state of the system, and record 
whether outcome "j" is obtained again, 

c) Repeat Steps a) and b) many times and denote the 
frequency of obtaining "j" in each Step by Uj and Qj 
respectively. 

(see below for a discussion of the number of repetitions 
required to estimate Uj and Qj to desired accuracies). 

Step 3: Repeat Steps 1 and 2 K times, where K is dic- 
tated by the desired accuracy and confidence in estimat- 
ing ub 

(see below for a discussion of the size of K) . 

Step 4: Compute the upper bound u6 of r^^^ (Eq. (7.14)) 
to accuracy and confidence dictated by the K trials 
{(li,l2),...,(ifi,if2)}fromStep3, 

ub^^^il-Z) 

- YTd V'l^T. V^^^VQfeiQfe. ■ (8-1) 

This concludes the protocol. 

We should emphasize various points about the proto- 
col. First, similar to the protocol with no output states 
(see Sec IV), the number of trials required to implement 
the above protocol is independent of d, and only depends 
on the desired accuracy and confidence of the estimates 
in each of Steps 2 c) and 3. See See's IX A and IX B for 
respective discussions about the time-complexity of each 
of these steps. 

Second, ub can be computed using only: 

1. Applications of the noisy measurement (sequen- 
tially on output states as well) and, 

2. The ability to prepare each of the d pure input 
states Ilj. 

Lastly, it is straightforward to show that the two neces- 
sary conditions given previously in Sec. V D for ub {lb) to 
be a useful upper (lower) bound also hold here. Indeed, 

1. In the limit of \, 0, 

ubiQ (8.2) 

and. 
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2. ub scales well in d. 

The first condition holds since as rg^M 4-0 it must be the 
case that Ek — > 11^ and pk — > Hk- Therefore, for every fc, 
Mfc — ^ 1 and Qk — > 1. The second condition holds using 
an analogous argument as that given in Sec. VD 2 



IX. RESOURCE ANALYSIS 

In this section we discuss the time-complexity and re- 
sources required for the protocols in detail, and also com- 
pare the time-complexity of the protocol with a full re- 
construction of the noisy measurement (ie. a full recon- 
struction of the POVM elements Ek). First, we analyze 
the number of trials required in Step 2b) of the proto- 
col in Sec. IV and Step 2c) of the protocol in Sec. VIII. 
Afterwards, we analyze the number of trials required in 
Step 3 of both protocols. 



A. Number of Trials Required in Step 2b) of 
Sec. IV and Step 2c) of Sec. VIII 

We first explicitly analyze the number of trials for Step 
2b) of Sec. IV. The discussion carries over in a straight- 
forward manner to Step 2c) of Sec. VIII. 



1. Step 2b) of Sec.IV 

For each k e {0, 1, d — 1}, we would like to under- 
stand how many samples are required to obtain an esti- 
mate of Uk- Let us fix k and define Zk = tr(nfe(l — Ek)) 
where, since {Ekjf^i is a POVM, 



(9.1) 



Then, 



Zk > 0, 
Uk + Zk =^ tr(nfe) 



1, 



(9.2) 



so {uk, Zk} forms a probability distribution on the binary 
measurement outcome space {"fc", "not fc"}. 

Let Uk be the estimator of Uk obtained in Sub-step 2b 
of the protocol. Since we are trying to find a parame- 
ter that defines an underlying probability distribution, 
we can use maximum likelihood estimation (MLE) tech- 
niques [51]. As the trials are independent and identically 
distributed, standard MLE gives the estimate of Uk to 
be equal to the frequency of obtaining outcome "fc", as 
in Step 4. Thus, if N is the number of trials and Uk is 
the number of times outcome "A:" is observed, standard 
MLE gives the obvious estimate. 



Uk 



nk 
N' 



(9.3) 



Since, in practice, Uk will be close to 1, one could uti- 
lize additive smoothing methods for MLE, which assign 
higher weight to low probability outcomes. In our case 
there is one low probability outcome "not fc". The for- 
mula for ilk under a A-additive smoothing method is. 



Uk 



nk + >^ 
iV + 2A' 



(9.4) 



where the factor of 2 multiplying A in the denomina- 
tor arises from there being two possible outcomes of the 
experiment, and is required to ensure {ufe, 1 — Uk} is a 
probability distribution. Note that when this smoothing 
technique is used, it will fairly bias the estimate of Uk to 
be smaller than the actual value of Uk, so the estimation 
will be honest. 

The key point in terms of time complexity is that the 
number of trials N required to estimate Uk to accuracy 
e with confidence 1 — (5 is independent of the size of the 
system. Since the set of possible probability distributions 
{ufc,Zfc} satisfies certain consistency conditions, the esti- 
mator Uk converges in distribution to Uk, 



. D 

Uk — > Uk. 



(9.5) 



Moreover, the mean of Uk is equal to Uk and the variance 
of Uk scales as 



Var {ilk) 



1 



NI{{uk,Zk}) 



(9.6) 



where I{{uk, Zk}) is the Fisher information [i'l] of one 
observation of the true probability distribution. The 
Fisher information contained in one observation is a mea- 
sure of how much information about Uk is gained on av- 
erage from observing { "fc" , "not fc" } (with distribution 
{ufc,Zfe}). If Ufe ~ 1 then the amount of information 
gained on average is large. 

Since {"fc", "not A:"} is a Bernoulli random variable, it 
is possible to explicitly compute the Fisher information 
of one observation, 



I{{Uk,Zk}) = 



1 



1 



UkZk Wfe(l-Ufc)' 

Hence, from Eq.'s (9.6) and (9.7), we have that 

Uk{l - Uk) 



Var (ufc) 



N 



(9.7) 



(9.8) 



Now, suppose we want to estimate Uk to accuracy e with 
confidence 1 — 6, ie. 

V{\uk-Uk\>e)<5. (9.9) 

By Chebyshev's theorem we have for any integer j > 1, 



'(|Ufe - Uk\ > j<j{uk)) < 



1 



(9.10) 
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Choose js to be the smallest j such that 

1 



< S. 



Then we have 



and so we set 



This gives 



js\/uk{l - Uk) 



N 



< e. 



N > 



j|ufe(l - Uk) 



In total, if TV > ^'^^""^/.'""^ with js > j then, 
V {\uk - Uk\ > e) < S. 



(9.11) 



I- I . -Uk)\ 1 , . 

\uk -Uk\> — 7= < (9.12) 



(9.13) 



(9.14) 



(9.15) 



Since N is independent of d, the number of trials re- 
quired to estimate Uk to accuracy e and confidence 1 — i5 
is independent of d. 



B. Number of Trials Required in Step 3 of 
Sec. (IV) and Step 3 of Sec. (VIII) 

Let us now discuss how many trials are required in Step 
3 of each of the protocols. The argument is the same 
for each protocol so, without loss of generality, we use 
the notation from Sec. IV (the discussion for Sec. (VIII) 
follows by replacing X with Z and making appropriate 
changes). We have that X is equal to the expectation 
value of the random variable X : 2? — > [0, 1] by 



X{l,m) = ^UiUr, 



(9.17) 



(for Sec. VIII this will be Z{l,ni) — y^uiUj^VQiQin)- 
Suppose one wants to estimate X to accuracy e and con- 
fidence 1 — 6 using K independent samples of X. Let 
Sk denote the random variable that represents the av- 
erage of K independent samples {Xi, Xj^} of X. By 
Hoeffding's inequality. 



¥{\S, 



K 



X\> e) < 2e'<<fi^ 

~2Ke? 



where [a, b] is the range of X (here [a, b] C [0, 1] 
setting 



(9.18) 
Hence, 

(9.19) 



2. Step 2c) of Sec. VIII 

In the protocol of Sec. VIII we have to estimate both Uk 
and Qk where, is the probability of obtaining outcome 
"fc" on the first measurement and Qk is the probabil- 
ity of obtaining outcome "fc" on a second measurement 
given the result of the first measurement is "/c" . Note 
that observing "fc" in the second measurement is also a 
Bernoulli random variable, so the discussion from above 
carries over in an analogous manner to estimating Q^. 

Let A^i and A'2 be the number of trials required to esti- 
mate Mfe and Qk to each of their respective accuracies and 
confidences. Suppose we perform the first measurement 
iVi times (so we have estimated Uk to its desired accuracy 
and confidence). If the number of times "fc" is observed 
over these trials is greater than N2 then we have esti- 
mated both Uk and Qk to their desired accuracies. If the 
number of times "/c" is observed in the first measurement 
is less than A'2 then we keep repeating until N2 outcomes 
of "fc" are recorded in the first measurement. Let M2 be 
the number of times the first measurement has to be per- 
formed before A''2 values of "fc" are recorded. The total 
number of trials is no more than 



gives. 



K 



ln(f)(fe-a)^ 
2e2 



< 



In I 



2e2 



(9.20) 



which is independent of d. In practice, 6 — a << 1 which 
will improve this bound on K . 



N = min{A^i, 1/2}- 



(9.16) 



C. Complete Reconstruction of POVM Elements 



For completeness, we provide an explicit protocol and 
determination of the time-complexity of complete tomog- 
raphy of the noisy POVM elements {Ek}k=i- The idea is 
to input various pure states into the noisy measurement 
and analyzing the frequency of obtaining particular out- 
comes. If the matrix representation of Ek in the |V'm) 
basis is written as E]:-' , then there are ^-5^ elements on 
the main diagonal and upper triangular section of Ek 
that need to be estimated (since each Ek > 0, the lower 
triangular part of Ek is completely determined by the 
upper triangular part). One can estimate the elements 
of Ek by first defining the pure states. 



As N only depends on the accuracies, confidences, and 
values of Uk and Qk, it is independent of d. 
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Then, since, 

(9.21) 

we have, 

~ (^) t^-(^fcn,) - tr(i?fen,). (9.22) 

The algorithm to determine the set {EkYl^i foUows. 

Step 1: For each state IV'jj') S {\'4'k,k)}t=n input IV'ij') 
into the noisy measurement £ many times and record the 
frequency of obtaining each of the d different possible 
outcomes "/c". This gives, for each fc, ir{Ek\il]j.j){'4'j,j\)- 

Step 2: For each state 1-0^^^) € {\i^k,i)}k,i=n input 
I'f/'^j) into the noisy measurement £ many times and 
record the frequency of obtaining each of the d dif- 
ferent possible outcomes "fc". This gives, for each fc, 
tr(i?,|0+.)(V'+|)=tr(i?fen+.). 

Step 3: For each state |V'i~j) € {l'0fc,i)}fc,i=ij input 
IV'Tj) into the noisy measurement £ many times and 
record the frequency of obtaining each of the d dif- 
ferent possible outcomes "fc" . This gives, for each fc, 

tr(£;,|7^,:^.)(V',^^.|)-tr {EuHT^). 

Step 4: Combine all of the elements estimated in Steps 
1 through 3 to construct the Ek ■ Step 1 gives the diagonal 
elements of the E^ since 

Ej:' ^ir{Ek\^,,,){i;,^,\). (9.23) 

Eq. (9.22) and Steps 1 through 3 give all of the off- 
diagonal elements i^^:"'. This concludes the protocol. 

The number of trials one will have to perform is again 
dictated by MLE. It is important to note that, as opposed 
to the simple Bernoulli procedure described in Sec. IX A 
for estimating the Uk^ the MLE procedure in this case is 
more involved. This is because when the noisy measure- 
ment is performed, one must keep track of which value of 
J is obtained (not just whether the outcome was j or not). 
Thus the number of trials in each step will be greater 
than that required to estimate the Uk because events with 
small probability may rarely be seen (if at all). Smooth- 
ing techniques will likely have to be employed to ensure 
rare events are not assigned zero probability. 

This discussion shows that the determination of each 
tv{Ek\(t>){(l)\), where |(/>) is one of IV'jj), or 
requires greater time-complexity than that of estimat- 
ing each Uk- Since there are d + 2^^fcii = such 10), 



one will have to estimate d^ different probabilities over 
Steps 1 through 3 using MLE. Hence, the full reconstruc- 
tion requires the estimation of d^ probabilities with more 
complicated post-processing of the measurement data (as 
well as a larger number of trials to estimate each prob- 
ability). In addition, there is an added complexity in 
preparing d^ different input states since more complex 
rotations may be required. 



X. DISCUSSION 

We have provided simple experimentally scalable pro- 
tocols for obtaining upper (lower) bounds on the average 
measurement error (fidelity) of quantum measurements. 
We discussed conditions for the validity of the bounds 
and explained why they should be valid in extremely gen- 
eral situations. The bounds could also be useful as actual 
estimates of the average measurement error (fidelity) in 
the small-error regime. We have presented a set of nu- 
merical examples for a single-qubit system. In every in- 
stance the bounds were shown to be valid (the bounds 
were not violated for any magnitude of coherence param- 
eters in the noisy POVM operators). In addition, as ex- 
pected, the bounds became better approximations as the 
magnitude of the coherence increased. This provides fur- 
ther evidence that the bounds should hold in extremely 
general situations and should be useful in practice. 

The protocols are scalable and only require the ability 
to prepare states from a basis set, and perform the noisy 
measurement (sequentially when there are output states) 
to estimate d probabilities. In addition, post-processing 
of the data is completely straightforward and avoids the 
difficulties in associating large sets of tomographic data 
to valid mathematical objects. This can be compared 
with a full reconstruction of the POVM elements of the 
noisy measurement, which requires the preparation of d^ 
input states, the estimation of d^ probabilities, and more 
involved post-processing. In addition, the d^ input states 
can be highly complex. In many situations, these input 
states will have to be prepared using complicated unitary 
rotations. 

There are a number of different questions and avenues 
for future research. First, we have focused on the case 
of rank-1 PVM's, however we expect that our results 
can be extended to higher-rank PVM's, especially low- 
rank PVM's in large Hilbert spaces. As well, since any 
POVM can be implemented via a PVM on an extended 
Hilbert space, our protocol can potentially give informa- 
tion regarding the quality of POVM measurements. It 
will also be interesting to analyze the extent to which the 
ideas presented here can be used to characterize non-ideal 
POVM measurements (as opposed to PVM's). While we 
have shown the bounds derived here should hold in large 
generality, a deeper understanding of the validity of the 
bounds will clearly be useful. 

Lastly, the algorithms given here require the ability 
to prepare the basis states {IfAfe)} which constitute the 
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measurement. In practice, these states typically have 
errors. We expect that running the protocol with noisy 
states will still provide valid bounds, however they will 
not be as tight as the upper bounds obtained when the 
protocol is run with perfect states. For instance, if the 
noisy input states are a convex combination of elements 
of {IV'fc)}) then the obtained bounds are still valid. 

The scalable protocols presented here can be useful for 
determining the quality of experimental quantum mea- 
surements. There is still much to investigate with re- 
gards to useful metrics for comparing measurements and 
proposing experimentally efficient methods for character- 



izing measurement devices. As experimental quantum 
systems scale to larger sizes, such methods will be indis- 
pensable for characterizing and controlling multi-qubit 
systems. 
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